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Abstract 

The temperature behaviour of meson condensates < ctq > and < erg > is 
calculated in the SU{3) x S'C/(3)-linear sigma model. The couplings of the 
Lagrangian are fitted to the physical vr, K, r], rj' masses, the pion decay con- 
stant and a 0^(1 = 0) scalar mass of nia = 1-5 GeV. The quartic terms 
of the mesonic interaction are converted to a quadratic term with the help 
of a Hubbard-Stratonovich transformation. Effective mass terms are gener- 
ated this way, which are treated self-consistently to leading order of a 
expansion. We calculate the light < qq > and strange < ss >-quark conden- 
sates using PCAC relations between the meson masses and condensates. For 
a cut-off value of 1.5 GeV we find a first-order chiral transition at a critical 
temperature Tc ~ 161 MeV. At this temperature the spontaneously broken 
subgroup SU{2) x SU{2) is restored. Entropy density, energy density and 
pressure are calculated for temperatures up to and slightly above the critical 
temperature. To our surprise we find some indications for a reduced contri- 
bution from strange mesons for T >Tc. 



* Supported by the Bundesministerium fiir Forschung und Technologie (BMFT) under contract 
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1 Introduction 

The study of finite temperature QCD is important both for theoretical and exper- 
imental hadron physics. Theoretically a hadron gas in thermal equilibrium is the 
simplest system to study the dynamics underlying deconfinement and chiral symme- 
try restoration. Experimentally a nearly equilibrated hadron gas with a transverse 
radius R± ~ 10 fm is supposed to give the dominant contribution at a later stage of 
relativistic heavy ion collisions. In principle this allows an experimental check of the 
equation of state of hadrons. In the following work we approach the phase transition 
region from the low temperature side. Below T ~ 100 MeV pions are known to be 
the most relevant degrees of freedom. Above this temperature region also heavier 
hadrons give a non-negligible contribution to the condensates and thermodynamic 
quantities [|l|. One way of including part of the heavier mesons is provided by the 
choice of SU (3) x SU (3) as chiral symmetry group rather than SU (2) x SU (2). The 
linear SU{3) x SU{3) sigma model includes a nonet of pseudoscalar (O^)-fields and 
a nonet of scalar (0+)-fields [g]. The spontaneous breaking of the SU{3) x SU{3) 
symmetry leads to massless (O^) Goldstone modes. Obviously a massless pseu- 
doscalar octet does not provide an adequate approximation to the experimentally 
observed meson spectrum. Therefore we include explicit symmetry breaking terms 
to account for the physical mass values of the octet-fields. A cubic term in the 
meson fields guarantees the correct mass splitting of the 77 — 77' masses which is due 
to the [/(l)-anomaly. It reflects the 't Hooft-determinant on the quark level. 

2 The model at zero temperature 

For a Euclidean metric the Lagrangian of the linear sigma-model is given as 
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^ = ia^$a^<l>+-i/i^tr$$+ + /,(tr<l>$+)' + /2tr ($$+)' 

+g (det $ + det + eoao + esc^s, (1) 

where the (3 x 3) matrix field is given in terms of Gell- Mann matrices \t [i = 
0, . . . , 8) as 

1 ^ 

* = ^E(^^ + ^^^)^^ (2) 
Here ag and vr^ denote the nonets of scalar and pseudoscalar mesons, respectively. 
As order parameters for the chiral transition we choose the meson condensates < 
ctq > and < erg >. The chiral symmetry of C is explicitly broken by the term 
{—Eoao — EsCTs), corresponding to the finite quark mass term rriqqq + m^ss on the 
quark level. The chiral limit is realized for vanishing external fields Eq and Es- Note 
that the action S = J d^xdrC with C of Eq. (|1]) may be regarded as an effective 
action for QCD, constructed in terms of an order parameter field $ for the chiral 
transition. It plays a similar role to Landau's free energy functional for a scalar 
order parameter field for investigating the phase structure of a $^-theory. First 
conjectures about the chiral phase transition were based on a renormalization group 
analysis in momentum space f^, ^. An e-expansion in d = 4 — e dimensions has 
been performed by lacobson and Amit 0] . When it is applied to the Lagrangian of 
Eq. (|l]) with g = 0,Eo = = Es, it predicts a first order chiral transition. For three 
fiavours the det-term is cubic in the field components. Hence a non-vanishing g will 
further support the first order nature of the transition. In contrast finite mass terms 
{eq ^ Es) may change the transition to a smooth crossover behaviour, if their 
values are large enough. 

Patkos and Frei confirmed the first order nature of the chiral transition in the 
chiral limit of the SU{3) x SU{3) linear sigma model [Q. The result was obtained 
in a saddle point approximation to the free energy functional in three dimensions 
(dropping the imaginary time-dependence of the fields $ in Eq. (|l])). In a subsequent 
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work IP it was shown that non- vanishing pseudoscalar meson masses actually change 
the first oder transition to a smooth crossover in the condensates, if otherwise the 
same approach is followed as in |p. 

In the present work we have extended the former approach to a treatment of the 
full four- dimensional theory, keeping all Matsubara frequencies in the effective po- 
tential. The reason is that a complete dimensional reduction from 4 to 3 dimensions 
may be approximately realized at high temperatures or for a second order phase 
transition. In the present work we are interested in the low temperature region; we 
also cannot expect a second order transition in the presence of explicit symmetry 
breaking terms. Therefore contributions from non-zero Matsubara frequencies may 
even qualitatively change the results. This is in fact what we will demonstrate in 
this paper. 

The six unknown couplings of the sigma-model (Eq. (||)) (/ig, /i, /2, g, Eq, Es) are 
assumed to be temperature independent and fitted to the pseudoscalar masses at zero 
temperature. Further experimental input parameters are the pion decay constant 
/tt = 93 MeV and a high lying (0+) scalar mass m„ = 1.59 GeV (cf. Table 1). 
For the remaining scalar masses and the coupling constants we obtain the values of 
Table 1. 



Input 


[MeV] 


rriKlMeV] 


m,j[MeV] 


TUr,' [MeV] 


U [MeV] 


rria- = fnfg [MeV] 


138.04 


495.66 


547.45 


957.75 


93 


1590 


Output 


fil [GeV^] 


fi 




9 [GeV] 


^0 [GeV3] 


^8 [GeV3] 


0.758 


12.166 


3.053 


1.527 


0.02656 


-0.03449 


rriao [MeV] 


m^. [MeV] 


mj/ [MeV] 


fK [MeV] 






914.05 


913.35 


764.71 


128.81 







Table 1: Tree level parametrization of the SU{3) x SU{3) linear sigma model (input 
data taken from experiment). 



The interpretation of the observed scalar mesons is controversial. There are good 
reasons to interpret the (O"*") mesons at 980 MeV as meson bound states [0. The 
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model underestimates the strange quark mass splitting in the scalar meson sector, 
the value for rriK* comes out too small. 

The effective theory can be related to the underlying QCD Lagrangian by com- 
paring the symmetry breaking terms in both Lagrangians and identifying terms 
with the same transformation behaviour under SU{3) x SU{3). Taking expectation 
values in these equations we obtain the following relations between the light quark 
condensates, strange quark condensates and meson condensates 



2m + nis Z[m — nis) 

<ss> = ^ [~''^ <ao> , <as>. (3) 

2m + mg [m — mg) 

We use m = (m„ + md)/2 = (11.25 ± 1.45) MeV and m, = (205 ± 50) MeV for the 
light and strange quark masses at a scale A = 1 GeV [If. From the scalar meson 
condensates at T = 0, (do) = 0.144 GeV and (ag) = -0.0415 GeV we get 



{qq) = - (235 ± 60MeV)^ 

(ss) = - (290 ± 30MeV)^ (4) 

in accordance with values from PCAC relations |Q within the error bars. Since we 
treat the coefficients Eq, Ss of < ctq > and < os >, and m, of < gg > and < ss > 
as temperature independent, we will use Eqs. (^) for all temperatures to translate 
our results for meson condensates into quark condensates. 

We also check that the pseudoscalar meson mass squares, in particular and 
m^- are linear functions of the symmetry breaking parameters eQ,es- Varying eo,es 
while keeping the other couplings fixed we can simulate the sigma model at unphys- 
ical meson masses. Since the current quark masses are assumed to depend linearly 
on aiid Es, an arbitrary meson mass set can be related to a mass point in the 
(m„ rf, ms)-plane by specifying the choice of (eoj^^s)- This may be useful in order to 
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compare our results for meson (and quark) condensates with lattice simulations of 
the chiral transition. 

3 Thermodynamics 

The thermodynamics of the linear sigma model is determined by the partition func- 
tion with the Lagrangian of Eq. (|I|) 

Z = JV(^expi^-J^^dTjd^xC{<^{x,T))Y (5) 

We will treat Z in a saddle point approximation. The saddle point approximation 
amounts to the leading order of a 1/A^-expansion in this model, where N = 2Nj = 
18. Note that C of Eq. (|^) would be 0{N)- invariant, if /2 = and g = 0. Our input 
parameters lead to non- vanishing values of /2 and g, therefore the 0(A^)-symmetry 
is only approximately realized. 

We calculate the effective potential as a constrained free energy density f/e//(^0) ^s), 
that is the free energy density of the system under the constraint that the average 
values of ctq and as take some prescribed values and ^s- The values ^Omin and 
^8min that minimize Uejf, give the physically relevant, temperature dependent vac- 
uum expectation values, i.e. < ctq >= C,Omin^ < ^8 >= ^8^i„- Hence we start with 
the background field ansatz 

o"o = ^0 + o-Q 

C^8 = ^8 + C^8> (6) 

where ctq and cTg denote the fluctuations around the classical background fields 
and .^8- All other field components are assumed to have zero vacuum expectation 
value, i.e. = cr^ for £ = 1, . . . , 7 and vr^ = vr^ for £ = 0, . . . , 8. The relation between 
the effective potential f/g// and Z is given by 
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. ^-l3VU,ff (^0,^8) 

J VM J ao{x, r) - ^0 5 J as{x, r) - 6 
crAx,T 



lis 
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(7) 



where ' /' is a short-hand notation for Jq dr J d^x. 



Next we insert the background field ansatz @ in C and expand the Lagrangian 
in powers of $' = Yle=o{'^e~^'^'^e)'^e- The constant terms in $' lead to the classical 
part of the effective potential Udass 



+ {fi + jj ^0 + (/i + d + 2 (/i + /2) eo^s - eo^o - Ss^s- 

Linear terms in $^ vanish for all ^ = 0, . . . , 8 due to the 5-constraints in Eq. (|^. 
Quadratic terms in define the isospin multiplet masses m^, where Q = 1, . . . ,8 
labels the multiplets. The contribution to the Lagrangian is denoted by C^'^^ 

^^'^ = lT.9iQ) {d,v'Qd,V>i + ^IfQf'i) ■ (9) 
^ Q 

Here ip'g denotes ctq for Q = 1, 4 and ttq for Q = 5, . . . , 8, g{Q) is the multiplicity 
of the isospin multiplet. We have g{l) = 3 for the pions, g{2) = 4 for the kaons, 
g{3) = 1 = g{A) for r],-)]', respectively. Correspondingly, the multiplicities for the 
scalar nonets are g{5) = 3, g{6) = 4, g{7) = 1, g{8) = 1 for the ao, -ft'o, /o, /q- 
mesons. Typical expressions for the masses are the pseudoscalar and scalar triplet 
masses 
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< = -^^l+{^f^ + lf2)eo+{^f^ + lf2)es 

+ lV2Mo^8 + -^g^o - 2^g^8 (10) 

< = + (4/1 + 4/2)^0' + (4/1 + 2/2) el 

+^Mo^s - -^g^o + 2\/|^?e8- (11) 

The cubic part in $' will be neglected, while the quartic term C^^^ 

£(4) = /i(tr$'<l>'t)2 ^ j2tr(<l>'<l>'t)2 (^2) 

is quadratized by introducing an auxiliary matrix field I](x, r). This is a matrix 
version of a Hubbard-Stratonovich transformation [^. We have the identity 



exp y ' dr J d^x{-) fi (tr ^ + /atr ($'$'1 

const ■ / VH exp — / dr I d x 

Jc-ioo ll6{e + 3a)'^ Jo 

{tr + 2fxltr S - 8(e + 3a) [eti 
+atr S ■ tr + //^£tr($'$'t) + 3a/x^tr ($'$'+ 



(13) 



where 



2£q; + 3a^ = fi 
^ h 



(14) 



In the saddle point approximation we drop / PS. As S'f/(3)- symmetric ansatz 
we use a constant diagonal matrix 



/ s 
S = s 
V s 

The choice of s will be optimized later, cf. Eq. 



(15) 



When Z is rewritten upon 
using Eq. (0), of Eq. (|l^) is replaced by given as 
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Hence the effect of the quadratization procedure is to induce an extra mass term 
(s + /ip) and a contribution t/gaddie to Ue/f, which is independent of and ^8- 

We are not aware of an analogous identity to Eq. ([T3|), which includes the det- 
term of C and leads to a tractable form. Therefore we drop the cubic term in $' as 
mentioned above. 

This way we finally end up with the following expression for Z 



g 

Q=l 

where 

X^Q = ml + fil + s. (19) 

Thus we are left with an effectively free field theory. The only remnant of the 
interaction appears in the effective mass squared Xq via the auxiliary field s. 

The choice of a self-consistent effective meson mass squared has been pursued 
already in Refs. |^. This is an essentially new ingredient compared to earlier 
calculations of the chiral transition in the linear sigma model . The positive con- 
tribution of s to the effective mass extends the temperature region, where imaginary 
parts in the effective potential can be avoided. In general, imaginary parts are en- 
countered, when the effective mass-arguments of logarithmic terms become negative. 
They are an artifact of the perturbative evaluation of the effective potential and of 
no physical significance, as long as the volume is infinite. In our application the 
optimized choice for s will increase as function of temperature and lead to positive 
Xq over a wide range of parameters. 
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Gaussian integration over the fluctuating fields $' in Eq. (0) gives 

Z{^0, ^s) = exp {-pV [[/class + f/saddle + 

4E«W)i:/(g5ln(,^(.;..J))]} (20) 

where 

cul^e + X^Q, (21) 

and 

ul ^ {2nn/Tf (22) 

denote the Matsubara frequencies. In contrast to our former approach p we keep 
all Matsubara frequencies and evaluate J2n£^ i'^ the standard way, see e.g. ||lT|. The 
result is 



^(^0,^8;^) = e-^^^-(«-««^^) (23) 

Uesi^O, S) = [/class + f/saddlc + [/th + f/o (24) 

Here we have indicated that Z and Ues still depend explicitly on the auxiliary 
field s. The integral in Eq. (^61) is regularized with a three-momentum cut-off A. 
The thermal contribution Utt vanishes at zero temperature and is finite for T > 0, 
while the zero point energy Uq diverges as A — >■ oo. 

The linear sigma model is a renormalizable theory, and the cut-off could be 
removed after a suitable renormalization prescription. Since we are dealing with an 
effective model, the need for such a renormalization may be less obvious. Anyway we 
do not believe in this model as an effective description for QCD, when the momenta 
exceed a certain scale, say A ^ 1 — 1.5 GeV. The necessity for a renormalization 
arises, when we postulate a matching between the physical masses and condensates 
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with the T = 0-values, and T approaches zero from above. Such a matching is 
guaranteed, if we impose the following subtractions on the zero point energy part 



Q^2 l™^hy. ■ ^^Q "^physJ + 2 d{m?f ^""^^y- ' "^physJ I- 

Here m^Yiys is given by trq of Eq. evaluated at =< > and =< >, 
i.e. for physical condensate values. The optimal choice s* for the auxiliary field s is 
then determined by 

"^^^^ ' 0, (28) 



ren 

s* 



ds 

where Ul^ equals U^^ of Eq. (0) with f/o replaced by U^'"'' of Eq. (^. 
Upon using Eq. (pSf ) it is easily verified that < >, defined as 

1 dlnZ 

is free of extra contributions from the zero point energies at T = if InZ = 
—(5VUIq^{^q,^^]s*). Thus a matching with < >r= is ensured. Similarly we 
find for the effective masses 

XI\t=q = + s + III = mjhys (30) 

for ^0 =< ^0 >, ^8 =< ^8 >, since s = — /Xq at T = 0. 



Note that the sensitivity to the cut-off in Eq. (|27|) is reduced from a A^-dependence 
to a 1/A^-dependence. This is a desirable feature in view of the uncertainties in a 
suitable choice for A. We have taken A = 1.5 GeV and kept the cut-off finite 
throughout the calculations. 

A further argument in favour of keeping the cut-off finite relies on results of 
Bardeen and Moshe |]T2[ on the 1/A^-expansion in 0(A^)-theories. According to these 



results a symmetry restored groundstate occurs as vacuum state at zero temperature, 
if the cut-off is sent to infinity. 
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Now we are prepared to determine the temperature dependence of the order 
parameters < > (T), < > (T) from the minima of f/efr'^(^o, ^s! s*). Thermody- 
namic quantities hke energy densities, entropy densities and pressure can be derived 
from Z in the standard way, if Z is approximated as 

4 Results 

For the parameters of Table |l| we vary the temperature and determine for each 
T the extremum of Ucs as a function of .^o; ^8 and s. The extremum is a min- 
imum with respect to and C,8 and a maximum with respect to s. (For zero 
temperature this is easily seen from the explicit form of f/ciass (Eq- ®) and t/gaddio 
(Eq. ([T7|)). Numerically the most convenient procedure is to determine the common 
zeroes in the derivatives and first, and then to check the desired 

minimum/maximum properties of the saddle point. 

In Figs. 1 and 2 we show the variations of ^-''^^\ and as a function 

° <qq>l =0 <ss>l =0 

of temperature obtained from < > (T) and < > (T) with the help of Eq. (^. 
We observe a gradual decrease of the light quark condensate, whereas the strange 
quark condensate stays remarkably constant. A first order transition occurs at 
Tc = 161 MeV for A = 1.5 GeV. Note that we have still a cubic term in the 
classical part of the potential (Eq. (|^)), the term we have dropped is cubic in the 
fluctuating fields. The critical temperature is determined in such a way that the 
pressure is continuous at T^. as a function of temperature. The values for the mesonic 
condensates show pronounced hysteresis effects, which are characteristic for a first 
order phase transition. 

At Tc the strange quark condensate does not drop to zero. Only the SU{2) x 
SU{2) part of the chiral symmetry is restored within numerical errors. While < 
(To >,< (Tg > have numerical errors A < a > / < a >= 0.1%, the error on 
the < qq > condensate is ^'l^^'^ = 4%. Systematic errors are attached to the 
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current quark masses. These errors give the main contribution to = 25%. 

An exact restoration cannot be expected, since the symmetry is exphcitly broken in 
the Lagrangian. When we derive the quark condensates from Eq. (^), we treat the 
quark masses and the external fields as temperature independent up to the transition 
region. This assumption may not be justified. The corresponding error is unknown. 

In our lowest order calculation of the effective potential we cannot distinguish 
between pole- and screening masses. By 'masses' we mean the effective masses 
Xq entering the arguments of the logarithm according to Eq. ([ToD , evaluated at 
the physical mass point =< ■Co > and ,^8 =< ■Cs >, cf. e.g. Eqs. ([lOD -(pTD. 
Thus the temperature dependence of these masses is determined by the temperature 
dependence of the condensates. The masses m^r and m/-/ become degenerate for 
temperatures T > within the errors ~ (139 ± 10 MeV), m^^ ^ (141 ± 10 
MeV) at Tc, when is approached from above. The degeneracy is a result of the 
vanishing light quark condensate. A zero value oi < qq> above Tc implies a relation 
between < ctq > and < erg >, which leads to ~ rn^,. Below Tc the meson masses 
stay remarkably constant. The pion mass and the scalar mass m^/ change from 

= (150 ± 10 MeV) and m^/ = (625 ± 10 MeV) at T = Tc approaching T^. from 
below to = (139 ± 10 MeV) and m^/ = (141 ± 10 MeV) at T = =161 MeV 
approaching from above. 

The behaviour of the scalar mass can induce an abrupt change in temperature 
dependent cross sections in contrast to a smooth variation in case of a second order 
transition or a crossover phenomenon. Temperature dependent cross sections may 
be realized in heavy ion collisions, when the transient quark gluon plasma cools 
down to the hadron phase and the hadron phase evolves until freeze-out. 

The /q above Tc can no longer decay into two pions, its width for T > Tc has 
a contribution from /q —>■ 27 decays with an invariant mass m?{2'-)) m^o- In 
the experiment one may see an anomalous amount of such (27) decays, when the 
system spends a sizeable time in the phase where chiral symmetry is restored. In 
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the experiment WA98 at CERN a measurement of the number of gammas to the 
number of charged mesons, i.e. mostly vr^, is planned. Above Tc, the n — K splitting 
is increased rather than reduced. The mass difference of Am = mK — mT^ = (400 ±20 
MeV) below (T < 161 MeV) is increased to Am = (516 ± 2 MeV) for T = 
161 MeV > Tc. Accordingly the strange meson contribution to the energy density 
in this temperature region is reduced compared to the low-temperature hadron gas. 

In Fig. 3 we give the energy density e/T'^ and pressure p/T^ as function of 
temperature for a cut-off A = 1.5 GeV. The gap in the energy densities at T^, which 
is a measure for the latent heat, is obviously rather small, about 10% of e at Tc. 
Sizeable contributions to e come mainly from 8 degrees of freedom, the pions, the 
kaons and the /o meson. The small value of the latent heat may be due to the 
vicinity of a (hypothetical) first order boundary in the (m^^, m7r)-mass diagram. At 
the location of this boundary the meson masses become so large that they wash out 
the chiral transition completely. We plan to check this explanation by investigating 
the chiral transition for unphysical values of the strange quark and light quark 
masses in the future. 

The existence of quasibound (^(/q) and vr modes may be a good approximation in 
the vicinity of Tc. Far above Tc the linear sigma model certainly fails as an effective 
model for QCD due to the lack of quark-gluon degrees of freedom. Nevertheless it 
would be interesting to study, at what temperature the full SU (3) x SU (3) symmetry 
is restored. At high temperatures the effective potential becomes proportional to 
X]qX^(Q)T^, the linear terms in the meson masses (Eqs. (p!0|)- ([Tl|) ) cancel and 
temperature tries to fully restore the broken symmetry. 

Finally we remark that our value for Tc is rather close to the lower limit of the 
Hagedorn temperature Th{Th ~ 160 MeV) |jl3[|. This may not be entirely accidental. 
In our 1/A^-expansion N means a large number of flavours, since 

N = 2-Nf. (32) 

In order to keep QCD an asymptotically free theory also the number of colours 
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Nc has to increase. Correspondingly our approximation is similar to Hagedorn's 
description of the hadron gas as a resonance gas. We expect that corrections from 
subleading terms in our 1/A^j-expansion will implicitly amount to corrections also 
to the large A^c-hmit. 
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Figure Captions 



(qq) 

Figure 1 Normalized light quark condensate 7— vs temperature. 

m)T=Q 

(ss) 

Figure 2 Normalized strange quark condensate — — vs temperature. 

Figure 3 Normalized energy density ^ and pressure ^ vs temperature. The 
decrease of these quantities above T fa 200 MeV as a function of T indicates 
the breakdown of our approximation scheme. 
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